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Abstract 



In this paper, an explanation of the Newton- Peiseux algorithm is given. 
This explanation is supplemented with well-worked and explained exam- 
I I pies of how to use the algorithm to find fractional power series expansions 

rn for all branches of a polynomial at the origin. 

< 

^ 1 Introduction 



Given a polynomial f{x,y), you can always solve for y in terms of x by means 
of a fractional power series 



> 

^^ The method for doing this is often explained in books in the context of a proof 

t^ of the algebraic closure of the field of fractional power series (e.g., [1]). Our pur- 

pose here is not to discuss the proof of such a theorem, but to explain clearly, 
with several examples, the algorithm for obtaining these expansions. 
Here is a general description of the algorithm. 



o 

00 

. . 1. Given f{x,y) — draw the Newton polygon of f{x,y). This is done by 

^ plotting a point on M^ for each term of / with the term kx°'y^ being 

k>< plotted to the point (6, a) The Newton polygon is the smallest convex 

;_] shape that contains all the points plotted. 

2. Take a segment of the Newton polygon from the set of segments where each 

point plotted is on, above, or to the right of the segments. 

3. The first exponent 71 will be the negative of the slope of this segment. 

4. Find /(x,x^i(ci-f 2/1)). 

5. Take the lowest terms in x alone. Since f{x,y) — 0, these must cancel and 

so you can solve for ci. 



6. Taking the values of 71 and ci and (3 — the 'x-intercept' on the Newton 
polygon of the segment we have chosen, we now find ,fi{x, j/i). 

8. Now repeat the process for fi{x,yi) to find 72 and 02- 

9. Continue this process until one of two things happen 

• fn{x, Hn) has a factor of y„. 

• The Newton Polygon of /n(a;,j/„) consists of a single segment with only 
two vertices, one on each axis. 

If the former occurs, either the y factor(s) can be factored out and you can 
continue, or (if you are not left with a segment in the Newton polygon) the 
series terminates. If the latter occurs, plug in the rest of the series, solving for 
coefficients by letting lowest powers of x cancel. Use the denominator of the 
slope of the segment to determine the increase in each 7^, (e.g., if the slope is | 
, each 7i, should be ^) . This will be greatly clarified by the following examples. 

Any polynomial of the form f{x,y) — can be solved for y in terms of x 
explicitly around a point. Since by simple translation any point on a curve can 
be moved to the origin, we will only expand a solution of f{x,y) = at the 
origin. The Newton-Puiseux algorithm gives us that any solution of f{x, y) — 
has the form y = cix^'^ + c^x^^^^'^ + c-^x^^^^'^^^'^ + • • • , where 7^ € Q+ and 
the Ci G Q"*". Any time 7^ is referred to, we are speaking of an exponent, and 
anytime a Ci is referred to, we are referring to a coefficient. 

2 Examples 

Example 1: In this example we show a more simple way to compute the ex- 
pansion that can be used if the second case in step 9, mentioned above, occurs 
where the Newton Polygon has a horizontal intercept of one. 

/(x, y) = 2x^ + x^y + Axy"^ + Ay^ = Q (1) 



We know that this polynomial has solutions of the form 

y = cix'^i + C2X^^+'^^ + C3x''i+''2+''^' + • • • . (2) 

The possible choices for 71 correspond to the negatives of the slopes of the lower 
segments of the Newton Polygon (Figure 1). The slopes of the lower segments 
of the Newton Polygon corresponding to (1) are -2, and -1. We will start by 
considering the case where 71 = ^ (^ 2) = 2. 

Let 7i = 2. Then by factoring x^^ (or x^) out of (2) we obtain, 

j/ = x^(ci +yi), (3) 



Figure 1: Newton Polygon for f(x,y) 




where ci is the first coefficient in the series and j/i is 



yi = C2X^^ + C3X 



p72+73 



C4X 



72+73+74 



(4) 



i.e. the rest of the series. 

So, substituting (3) into (1) we get, 

2x'^ + x4(ci + yi) + 4a;5(ci + y^f + 4x^{ci + y^f 







(5) 



The vertical intercept of the segment we used to find /i on the Newton 
Polygon, which we call /?, gives the terms of lowest degree in x alone. These 
terms must sum to zero because /(x, y) = 0. In this case, /3 = 4, so 
2x'^ + cix^ ~ and solving for ci we obtain ci = —2. Substituting ci into (5) 
and dividing by x* gives, 

x-^f{x, x2(ci + yi)) = 2 + (-2 + j/i) + 4a;(-2 + yif + 4x2(-2 + y^f 



2„,3 



J/1 + 16x - 16xyi + Axyt - 32x^ + 48x^yi - 2Ax^y^ + Ax'^yf 



(6) 



Figure 2: Newton Polygon for /i(x,j/i) 



We call this new polynomial /i(x, yi), so now that we have our new function 
fi{x,yi) our next goal is to find C2. We find the new Newton Polygon for 
fi{x,yi) (Figure 2), from which we get /32 = 1 and 72 = 1- Again, we let the 
terms of lowest degree in x alone equal zero. We have now reached a point in 
our algorithm where the calculations become much simpler. Notice that in the 
previous Newton polygon the intercept with the horizontal axis was one and we 
had only one segment. This tells us that we can calculate the rest of the Cj and 
7i directly from the current function, in this case /i.This is true because when 
we obtain such a Newton Polygon, we know that the powers of x are going to 
increase by at least 1/ (denominator of 72) from now on (with the possibility 
that Ci could be = ); in other words, 73 = 74 = • • • = 7n = • • • =. Since 72, 
we know that 



yi 



C2X + C^X^ + C4X^ + C^X^ + 



So, substituting yi into (6) and letting the lowest degree terms of x cancel we 
see 

C2X + 16x — 
so, C2 = — 16 

We find C3 in the same manner. To do this, we let the x"^ terms equal zero. 
Substituting yi into (6) we see that, 

C3a;2 - 16c2x2 - 32x2 ^ q 

C3 + 256 - 32 = 

so, C3 == -224 

Now, we can let the x'^ terms in /i be equal to zero. So, 

ax^ - 16c3X^ + Ac^x'^ + 48c2X^ = 

a + 3584 + 1024 - 768 = 

so, Ci = -3840 



Next, we let the x** terms equal zero. 

cgx" - 16c4X^ + 8c2C3X^ + iSc^x^^ - 2\c\x^ = 

C5 + 61440 + 28672 - 11184 - 10752 = 

so, C5 = -39504 

This process can be continued to calculate more terms, but we will stop with 
C5. So, one explicit solution of the polynomial is 

y = -2x2 _ ;^g^3 _ 224x4 _ 3840x5 - 39504x6 _,_... 

We have examined the case where 71 = 2, and now we must return to the 
original polynomial /(x,y) and examine the case where 71 = 1. Recall, 
/(x, y) ~ 2x4 _|_ ^2y _|_ 42,y2 _^ ^y3 _ Q^ L^^ -^^ — 1. Then, y — x{c\ + y\) and 
substituting 71 and y into (1) we obtain 

As in the case where 71 = 2, 

yi =: C2X''2 + C3X''2+''3 _^ c4X''2+''3+''* H 

From the Newton Polygon (Figure 1) we see that ji\ — 'i. Therefore, we can let 
the x^ terms equal zero. 

cix^ + 4cf x^ + 4ciX^ = 

ci(l+4ci+4c2) = 

ci(l + 2ci)2==0 

C\ U, — 2 , ~ 2 

Ci = is trivial, so assume, ci = — |. 
Substituting ci into (7) and dividing by x^ gives: 

x~^j{x, x(ci + yi)) = 2x + (-- + ?yi) + 4(-- + yxf + 4(- ^ + Vif 
= 2x + yi - 4yi + 4y2 + +3^^ - 6y^ + A.yl 
fi{x,y)^2x-2yl+4yl 

72 = 2, so 

yi = X2 (c2 + 2/2) so, 2x - 2x(c2 + j/2)^ + 4x3 (C2 + 7/2)^ = 

where C2 is the second coefficient and 2/2 is the rest of the series. Substituting 
yi into Ji{x,y) yields 

2x - 2X(C2 +y2f+ 4x1 (C2 + j/2)^ = 

From the new Newton Polygon (Figure 1) we see that /32 = 1, so we can let 
the X terms equal zero. 



2x - 2xcl = 

C2 = ±1 (8) 

Let C2 = 1. (Wc must return later and consider the case where C2 = — 1-) Then 
substituting C2 into fi{x, y) and dividing by x we get 

x- V(x, y)i = 2 - 2(1 + 2/2)' + 4x5 (1 + y^f 

= -4y2 - 2?/2 + 4a::5 + 12xUj2 + I2x^yl + 4x^2yl 
Where wg = cax'^^ + Cixi'-'+'>* + c^x'^'-'+f^+'i'' -\ 

1 3 

The horizontal intercept is one, so j/2 = c^x^ + C4X + C5X2 + • • • 
Letting the lowest degree terms of a; (xs ) equal zero we find 

-4c3 + 4 = 

Now, we can let the x terms equal zero. 

-ic^x - 2c\x + I2C3.T = 

-4c4 -'2 + 12 = 

so, C4 = I 

We now let the x^ terms equal zero. 

-4c5x5 - 4c3C4X5 + I2C4I + I2c\xi = 

-4c5 - 10 + 30 + 12 = 

so, C5 = 8 

As in the case where 71 = 2 we could continue this process and calculate more 
terms. Stopping here, another exphcit solution of /(x, y) is 

y=-5X + x5+x^ + |x5 + ^x^ H 



Now, we must return to (8) and consider the case where ci = —1. Substi- 
tuting ci into /i(x,y) and dividing by x'^ gives 

x-^h{x, y) = 2 - 2(-l + y2f + 4x5(-l + y^f 

= 4j/2 - 2^2 - 4x3 -I- 12x5y2 - ^'^x^vl + '^x^vl 

Let y2 = cgx'^--' + C4x''3+'r4 + c^x'^^+'^i+is _^ ^ ^3 = 1^ go 2/2 = csx^ + C4X + 

c| + • • • . Letting the x 2 terms equal zero, 

4c3X2 + 4x2 — 
so, C3 = -1 

Letting the x terms equal zero, 



4C4X — c§x + I2C3X = 
4c4 - 1 - 12 = 

so, C4 = ^ 

Now we let the x^ terms equal zero and solve for C5, 

4c5 + 12c4 - I2cl == 
4c5 + 39 - 12 = 

27 
SO, C5 = -^ 

So, we find that the last explicit solution for f(x,y) == is 

y = —\x - xi — x"^ + ^xi — ^x^ 
So, the three solutions for /(x, j/) = are 

y = -2x2 _ iQ^3 _ 224a;4 _ ^g^Q^^ _ 39504x6 

y = -^x + x^ + x^ + jx^ + ^x^ -\ 

y = —\x - xi — x^ + ^x5 - ^x^ 



Example 2: In this example we see that the Puiseux algorithm can be used to 
compute highly accurate fractional powers and expansions with imaginary and 
complex terms. 

Figure 3: Newton Polygon for /(x,y) 




/(x, y) = x'"^ + Sx'^ - 2x2y2 _ y3 _^ ^y'^ = 
We know that this polynomial has solutions of the form 

y := cix''^ + C2X^^+'''' + C:iX^^+^''+'<'' H . 



(1) 



(2) 



7i corresponds to the negatives of the slope of the lower segment of the Newton 
Polygon (Figure 3). The slope of the lower segment of the Newton Polygon 
corresponding to (1) is — |. So, 71 = — (— |) = |- 

Let 71 = I Then by factoring x"'^ (or a; 3) out of (2) we obtain, 

2/ = a;t(ci + yi), (3) 

where ci is the first coefficient in the series and 

yi = C^x'^^ + C3x''^+''3 _^ ^^3,72+73+74 _^ . . . (4) 

i.e. the rest of the series. 

So, substituting (3) into (1) we get, 

x^ + Sx-* - 2a;¥ (ci + y^f - x'^{ci + y^f + 2a;^ (ci + y^f ^ (5) 

The vertical intercept of the segment we use to find 71 on the Newton Poly- 
gon, which we call (3, gives the terms of lowest degree in x alone. These terms 
must sum to zero because /(x, y) = 0. In this case, /3i = 4 so, 8x^ — c\x^ = 
and solving for ci we obtain 

ci ==2,-l + i\/3,-l-i\/3 (6) 

We assume ci — 2, but we must consider the other two cases later. Substituting 
Ci into (5) and dividing by x^ gives, 

x"''/(x, x^{ci + yi)) == /i(x, yi) 

= X- 8x3 - Sxsyi - 2x5y^ - 12yi - Qyl - yf + 32x3 
+ 64:xiyi+4:8xiyf + 16x^yl + 2x^y* (7) 

We call this new polynomial /i(x, yi). So, now that we have our new function 
/i(x,yi) our new goal is to find C2. We find the new Newton Polygon (Figure 
4), from which we get /32 = | and 72 = |. 

Similar to (3), we factor x'^^ out of (2) and obtain 

yi=x3(c2 + y2) (8) 

where C2 is the second coefficient in the series and 

Subsituting (8) into /i we see that 

O^X- 8x3 _ 8x5 (c2 + J/2) - 2x2(c2 + y2f - 12x5 (C2 + y^) - 6x^3(02 + y2f 

- x^{c2 + y2f + 32x5 + 64x2(c2 + 2/2) + 48x5 (c2 + y2f + 16x^5^ (c2 + y2f 
+ 2x\c2 + y2f (9) 



Figure 4: Newton Polygon for /i 



• • • • 



bCJL 



012345 



P2 gives us the terms of lowest degree in x alone, which we can again set equal 
to zero. 

-8a;3 - 12c2a;3 = 

so, C2 = -| 

Substituting C2 into (9) and dividing by X3 we see 



_2 2 1 32 10 1168 4 

X 3/2(x,a;3(c2 + y2)) ==a;3 + — a;3 - __x3 

196 



128 



64 



104 



-X3 



-X3y2 - I2y2 -6x3y^-x3y-i -&Ax^y2 

99 64 8 89 80 64 10 

48a;'y2 + -jx^'y2 - 32x5y| + 16x3y^ -^ —x ^ y2 



16 1 



^„.2 



16 



^„,3 



II) 4 



y 2; 3 y^ - —X 3 y^ + 2x 3 y* 



(10) 



Now, we can see from Figure 5 that 73 = = , so 

2/2 = a;s(c3 + y3) 
where C3 is the third coefficient in the series and 



(11) 



Figure 5: Newton Polygon for /2 






We substitute y2 into /2 and find 

1 32 10 1168 4 128 8 64 , 104 a 196 5 , 

=.x3 H X 3 a;3 x^ -\ x -\ x-^ H x'^ (C3 + %) 

81 27 27 3 3 3 '^ ^^^^1 

- 12a;3 (c3 + 2/3) - 6a;5 [cy, + ?y3)2 - a;? (03 + 7/3)3 - 64a;5 (03 + y^) 

+ 48x3 (C3 + 2/3)2 _^ — x3(c3 + y^) - 32a;^ (C3 + ya)^ + 16x^(23 + y^f 

64 11, , 16/, ,9 I613, ,Q „ii/ na 

- ^x 3 (c3 + 2/3) + ya;*(c3 + ya) - -jx 3 (03 + ya)-^ + 2x 3 (ca + ya) 



(12) 



Since /2 = the X3 terms must also equal zero, so we can let 

a;3 — 12c3X3 = 
so, C3 = Y2 

Substituting C3 into (12) and dividing by xs we see that 



X 3/3(a;,x3(c3 + y3)) 



61 
324' 



-X 3 



119 



-a;3 



9353 

' 216 ' 
49 



-a;3 



1 

X 

324 
104 

196 



4 9217 2 

— a; 

1728 
x5 - 122/3 - 



14 

\ 5 

81 

64 



1 

27 
1 



X 3 



16 



-X3 



-X3 



10368 



-X 3 



64 



16 Q I4 44,^„„, ^ iv— - 

yx-'ya- g^^ ^3 - ^a;''?/! + —x^ya + 8x3^3 +48xSy^ + yx^ya 



55 10 in 9 8 ii 16 n 9 

— X 3 y3 + 4x 3 y^ + -a; 3 y3 + —X 3 2/^ 



1 

216 



X 3 y3 



1 
12 



13 9 



+ 3^; 3 y| - xya - 6xy| 



3073 



1 



48 



-X ya 



„2„,3 



3„,2 



a^'ya - a;^y^ - 32x3y^ 



16 



+ 16x^y3^-— x*y3^ + 2x-y^ 



(13) 



10 



Figure 6: Newton Polygon for /a 



Now, we can see from Figure 6 that 74 = | so we have 

ya = x5(c4 + y4) 

where C4 is the next coefficient in the series and j/4 is the rest of the series. 
Substituting j/3 into /a gives 

„ 2 16 o o, ,0 64 o, , 3073 7 

^-x 3 (c4 + y4r - x-\Ci + y^Y + —x^a + 7/4) - ^^.t3(c4 + Va) 



— 32x 3 (c4 + j/4)^ — 6.T ^^^ , y^j , 

o 



5 / n9 16 13 9 4 14 / 

3 (C4 + 2/4)-' + —X 3 (C4 + ^4)^ - -X 3 (C4 + i;4 

- -x^ (c4 + ^4) + 2x^ (c4 + ^4)'' - -xt (C4 + y4)2 - 12a;i3(c4 + 2/4) 



9353 



+ 16a;T(c4 + y4)^ 



+ 48x3(c4 + y4)^ + 4x''(c4 + y4) 

216 

16 10 196 5 1 11/ N 

yx 3 (c4 + y^) + — a;3 (C4 + y^) + — x 3 (C4 + y^) 



+ 8x3(c4 + y4) 



16 5, 
- —X (C4 + y4 



-xt(c4 + y4) ^ 
119 



5/ \2 ^^ — / \ ^ 4/ \ 

—X (C4 + yi) - —X 3 (C4 + 7/4) + -X (C4 + ^4) 



1 4 9217 o 14 . 1 11 16 8 61 10 

X X H X -\ X 3 -I x3 X 3 

1728 81 27 9 324 



1 



7 49 4 104 1 64 5 

X3 -\ X3 -\ X3 -\ X3 -\ X 

9 3 3 10368 



(14) 



From the Newton Polygon (Figure 6), /34 = |, so we let the X3 = and solve 
for C4. 



i^X5-12c4X3 =0 
SO, C4 = -^ 



11 



Substituting C4 into (14) and dividing by 

-ij-/ i/ , NN 23 4l32 

40901 ^2 



X3 



40119049 
93312 ^ ^ 216 



5047 10 
X ^ Va 



-X 



5053 o 



27 



35152 r 

-x^ 



319510 8 

x^ 



66317 10 

T 3 -|- 



729 '^'' 243 
1352 . 16 , 



324 

1847 



144x4y| 



■X3 



^ 81 
-2:2/4 - 



11645 4 
-X 2/4 



OU I O 9 zl '^ 

^^x 2/4 + 16a; y^ 
2704 



279695 M _ ^ 

2187 27 ^* 

1663 13 

X 3 W4 — 12t/4 

216 ^ ^ 

913952 16 118 

- X 3 4- 

6561 3 

, 11 1 s. 1. 
-4x3 yl-x^yl 



•4-3_6xty2 
140608 



+ 160x5 

+ if'vl 

1 I 2 
4^' 2/4 



2187 
11647 11 

+ ^2r^^ 

9353 2 52 . o 



16 14 o 

yx3,,- 

10 o 

-32x^2/4 



- 2x 3 y| 



729 



X 3 2/4 



9217 5 3601 13 

2; 3 _ X 2 

1728 324 



a;%4 



7388 



X3J/4 



- 24x 3 j/4 

59 7 

y^'2^4H 

10777 14 

X 3 y. - 

81 '^ 



92 4 

- y-^J/4 
208 16 o 

T 3 l/'^ 



9 
1661 

1^ 



x'yt 



14 9 

a; -'2/2 



(15) 



Figure 7: Newton Polygon for /4 

6 



tx. 



3 4 5 



It can be seen from the new Newton Polygon (Figure 7) that 75 = |, so 

2/4 = a;s(c5 + 2/5) 



12 



Substituting ^4 into /4 wc find that 

„ 23 ., ,2 92 2. ^ 35152 , 40119049 4 

=y X + 4.'^(C5 + ,5f + -x^(c, + ,,) + ^^x^ + ^^^^^^x^ 

40901 2 5053 , 2, , 11647 11 319510 g 

66317 10 1847 7 ,^^4 9353 2 279695 i4 

X 3 H x3 -|- 160x3 x3 X 3 

243 81 216 2187 

+ 144x 3 (c5 + ys)-' - 6x5 (c5 + y^Y - X 3 (c5 + y^f - — r^aj-^Cs + ys) 

216 

16 11, , 111, n9„.13/ s 

- -^x 3 (c5 + ys) - ^a; 3 (cg + ys)^ + 24x 3 (cg + yg) 

10777 16, ^ 2 7, ,0 1352 17, 
^^x 3 (c5 + ys) + gS; (cg + ys) + -^2; 3 (cg + y^) 

140608 fi, , „„ 14, ,0 3073 «, 

+ ^^Q x^icb + ys) - 32x 3 (eg + yg)^ - -^x^' (eg + yg) 

59 q, , 52 19 9 208 22 q s 

+ ya;^(c5 + ys) + ya; 3 (cg + yg)^ + — x 3 (cg + yg)"" - x5(cg + yg) 

1 1 Q 1 (-J(-J1 

+ 2x«(c5 + yg)^ + ^-x4(cg + yg)2 - — — x6(c5 + yg)^ + 16x6(cg + yg)^ 
3 36 

2704 20, ,2 11645 14, , 5047 ,, 

+ -1^^ ■' (C5 + ys) + ^^^ ' (cs + ys) - -i^x (cg + yg) 

7388 10, , 16 20, ,q 4 17, ,2 9217 5 

+ -i^x 3 (eg + yg) - yx 3 (eg + yg)-^ - -x 3 (cg + yg)'' - -^^^x^^ 

3601 13 913952 w 

We can see from the Newton Polygon (Figure 7) that /3g = |, so we can set the 
X3 terms equal zero and solve for cg. 

^xi - 12cgxt = 

^*^' ^^ ~ 2592 

This process could be continued to calculate more terms, but we will stop here. 
So, an explicit solution for /(x, y) is 

9353 10 26 8 I72, 4 

^=2592"-^+y"' + T2"' + r +'"' 

We must now examine the case where the solution to ci is imaginary. There 
are two imaginary solutions, but we will only show the process for one of these 
cases because they are quite similar. Let ci = — 1 + ^/2i. We substitute ci into 



13 



(5) and divide by x^. 

x~'^fi{x,x'^{ci +yi)) = .fi{x,yi) 



iiVSxi - SiVSyl - Idx^ - 2ix^ 



vl 



2x 3 j/i + 64x 3 j/i - 24ix 3 3 ^ y^ _|_ Qi^/^y^ - Ai\/'ix 3 yi 



4x3yi - 2x3yJ + 16z%/3a;3 
4x3 



8i%/3x5y3 



2/? 



8xty3 
4r 
6yi 



3y? 
(17) 



Figure 8: Newton Polygon for /i 



lO 









We can see from the new Newton Polygon (Figure 8) that 72 = |, so 

yi = xi{c2 + y2) 
where C2 is the second coefficient in the series and 

2/2 = C3x''3 + c4a;''3+74 _^ ^^^i^+n+ia _, 

We substitute j/i into /i to get 

0=a; + 4?V3a;3 - Ai^/?,x^ {c2 + y2) -16x3 - 24xi(c2 + ^2)^ - 8x^(c2 + 2/2)^ 

+ 2x*(c2 + y2f + 64x2(c2 + ^2) + QiVixi{c2 + ya) - SiVSxt (c2 + 2/2)^ 

- 24i\/3xi(c2 + y2f + 4x5 (c2 + j/2) - 2x^(02 + y2)^ + 16i\/3x3 

+ 8iV3x^ (C2 + 2/2)^ + 6x5 (c2 + 2/2) + 3xt (c2 + ^2)^ - x^{c2 + 2/2)^ + 4x3 

(18) 



14 



From the Newton Polygon (Figure 8) , /32 = | , so we can let the x 3 terms equal 
zero and solve for 02- 

6c2X3 + 6i^/3c2X3 + 4:i\/3x3 +4x3 — 
2 



Now, substituting C2 
x^\f2{x,xri{c2 + y2)) 



so, C2 

into (18) and dividing by X3 gives 



32 10 1168 4 

— rr 3 2; 3 

27 



- a;3 



-62/2- 



- 32a;2y2 



32 , 64 8 32 

64 10 16 10 2 

7^X 3 y2 + —X 3 ^2 



„ . 9 9 16 10 o 64 10 

24x^2/9 ^a; 3 y^ - -—x 3 1/2 -, „ 

y2 3 y2 27 *'^ 3 

32 8 ^ „ 3 9 196 4 „2o 4o ^S'^ 

- ya;5y2 + 16a;3y^ + — x3y2 + 3a;3 j,| -x^^y^- 8x^y^ 

'-^ 242\/3.T^y2 + 32iV3x'^y2 + —iV3x^'y2 

o 

„l„,2 , 52. /-^2 64 /- 8 32 /- 2 



-2x 3%^ - 3iV3x3j/^ 



52 2 

T 3 

3 



- SiVSx^yl - WiVSx^yl + —iVSx^ - —iVSx^ - — zVSx^ 

- 6i\/3j/2 



(19) 



Figure 9: Newton Polygon for /2 



• • • 






012345 



From Figure 9 we see that 73 = 1 1 so 

2/2 = a;s(c3 + j/3) 
where C3 is the next coefficient in the series and 2/3 is the rest of the series. 
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Substituting y2 into /2 gives 
„ 16 4, ,2 32 2 64 8 32 10 1168 4 i 

=y^ ("^a + y^y ^ y^ + ^^' + sT"' ' " ^r"^' ^ "^^ 

52 2 1 7 o 64 11 32 o 

- yX3 + 6X3 (C3 + yg) _ .^3 (c3 + y^Y - — x 3 (C3 + yg) - y.^ (c3 + Vz) 

- 24x1 (C3 + y:if + 2x^(03 + ^3)4 _ gxir (03 + y^f + — a;s (cg + ^3) 

+ 16a;^ (c3 + ^3)2 + 32x5 (C3 + ^3) + 3xt (C3 + ^3)^ - y x^ (03 + y3)3 

^ 1 , , 52 /- 2 64 ^ 8 32 ^ , 
+ 6iv3x5(c3 + ys) H iv3x5 JV3x5 i^Jix^ 

00 

+ yi\/3x3(c3 + yg) + 8i\/3x^(c3 + y3.f - 16^^3x^(03 + ^3)^ 

- 24i\/3x5(c3 + y^f - 3i\/3xt(c3 + y^f + 32^^3x5(03 + y^) (20) 



We can see from the Newton Polygon (Figure 9) that (3^ — ^, so we can let the 
X3 terms equal zero and solve for C3. 

6C3X3 + 1x3 + 6i\/3c3X3 = 



16 



Substituting C3 into (20) and dividing by a; 3 gives us 



1/3(0;, X3 (03 + ^3)) = 



USiVSx^yl SiVS: 



128x3 



■a^'yl 



64i\/3a;3 



2/3 



(*\/3- 
244xTr 



1)4 



(iy3 + l)4 

4 

1568x3^3 
3(i\/3 + l)4 

3073x2^3 

6(^^/3 + 1)4 

8x^y| 



9(^^/3+l)4 
t 2 



128a;TT'2 



81(i\/3 + l)4 
128a; Try3 



3(iV3+l)4 
96^3 



9(i\/3- 
8a;j/4 



1)4 



(2\/3+l)4 

1904x5 
27(i\/3 + l)4 
96xy| 

(iW+I)^ 

4x2 y2 



(i\/3+l)4 (zV3+l)4 («y3 + l)4 
32/3 8x~ 



7 

64x3 1 



8iy3x4y3 



:\/3 + l)4 (^^3 + 1)4 27(i\/3 + l)4 9(i\/3+l)4 

2x4 784^1 



(za/3 + 1)4 
128iy3x4y| 



13 
X 3 



3(^73+1)4 

512x3y| 
(n7f+I)4^ 

64x4 j/| 
3(i\/3 + l)4 



X 3 ZX ('84X3 

648(^^3+1)4 81(i\/3 + l)4 9(^73+1)4 
440x1^7/, 384xs,;2 128x3^3 

3(i\/3+l)4 



^J/3 



8x4^ 



440x 
27(^3+1)4 (j\/3 + l)4 
128x4y| 
3(zV3+l)4 

13 



32x¥ 



yi 



256x 



10 o 

3 2/3 



9(i\/3 + l)4 27(i\/3 + l)4 



13 9 

_ 2X3 y^ 

3(1^3 + 1)4 
512iV3xi 

3(i\/3 + l)4 

_ 96iV3y3 
(jy3 + l)4 
3073i\/3x2y3 

" 6(iV3+l)4 
128^^3x3^3 

~ 3(iy3 + l)4 
1664x3 
3(i\/3 + l)4 



2iV3x^ 
81(i\/3 + l)4 

9353ix\/3 
27(i\/3 + l)4 
2i\/3xTrj/| 



(l\/3+l)4 (2^/3+1)4 

16x"3" 
i\/3+l)4 

128i\/3xi 
9(n/3nj^ 

112i\/3x3 



32x " y| 
3(i\/3 + l)4 
9217i\/3x2 
216(i\/3+l)4 
8z\/3xT 



81 (V3- 

/3xTry3 



1)4 



27(i\/3+l)4 

^.tyo.^ ■■ 1^3 _ i-yo^ - y3 _ 16iA/3x^y| 

(iVS + 1)4 27(i\/3 + 1)4 {iy/3 + 1)4 
1568i\/3x3 UOiVSx^ 



3(i\/3+l)4 

112x3 

81(i\/3+l)4 

9353x 
27(i\/3 + l)4 



27(i\/3+l)4 

1024xtj/3 ^ 
3(«\/3 + l)4 

512xt 
3(^^/3+l)4 



2/3 



32i^/3x 3 
^ (z\/3 + l)4 
9217x2 
216(i\/3 + l)4 
384^ 73x5 7/3 

(iA/3+l)4 ^ ^ 
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By continuing this process it can be found that C4 = — — . ^ and 



_ 9353 

^ 1296(l+i\/3) 



So, a second exphcit sohition for f{x,y) is 



/- 4 29 X3 52a:; 3 9353x 3 

y = —1 + zv 3a;^ X 1 1 

3' 6{iV3 + l) 9{-l + iV^) 1296(i\/3+l) 

Using a similar process, the exphcit solution of f(x, y) in the case where 

ci = —1 — iyi can be found. So, we have that the explicit solutions for /(x, y) 

are 

9353 10 26 8 1 



2592 9 



X3 H X3 + -X +2x3 + . . . 



12 



-1 + i^/3; 



a;3 — -x 



52xt 



9353x1? 



^N 4 2 9 x3 

j/= (-1 -iV3)x3 - -X + 



3 (6(i\/3 + l) 9(-l + iV3) 1296(i\/3+l) 

52xt 9353x^03 

3^"" 



(6(-l + i\/3) 9(i\/3+l) 1296(-l + i\/3) 



Example 3: 



Figure 10: Newton Polygon for /(x, y) 

6 




In this example, we find that multiple terms of a Puiseux expansion can be 
calculated before the Puiseux jets differ. Here, the jets will appear the same for 
the first three terms and then split at ^. 

/(x, y)^y^ + 2x^y + x^ + x^y^ + xy^ + -/ + x^y + x^y^ - -x/ - -y^ (1) 
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We know that this polynomial has solutions of the form 

y = cix'<^ + C2x'^^+''^ + C3a;''i+''2+''^' + • • • . (2) 



The possible choices for 71 correspond to the negatives of the slopes of the 
lower segments of the Newton Polygon (Figure 10). The slopes of the lower 
segments of the Newton Polygon corresponding to (1) are both -2. We will start 
by considering this case where 71 = ^ (^ 2) ~ 2. 

Let 71 = 2 Then by factoring x^^ (or x"^) out of (2) we obtain, 

y = x^{ci + yi) (3) 



where ci is the first coefficient and yi is 

yi ^ C2X'^-^ + C3x''^+''3 + CiX'^^+'^^+'^* + • • • (4) 



i.e. the rest of the series. 

So, substituting (3) into / we get, 

=x''(ci + yif + 2x''(ci + yi) + x"^ + x^{ci + yif + x\ci + yif 

+ Jx«(ci + yif + x\c^ + yi) + x\c^ + y^f - ]^x\c^ + y^f 

-\^'\c,+y,f (5) 



The vertical intercept on the Newton Polygon (Figure 10), which we call 
/3, gives the terms of lowest degree in x alone. These terms must sum to zero 
because /(x, y) ~ 0. In this case, /3 == 4 so, x^c\ + x'^ci -\- x'^ — Q. The x" part 
of the terms can be divided out and will not be included in subsequent steps. 
Solving for ci we obtain ci — —I Substituting ci into (5) and dividing by x* 
gives, 

x^^f{x,x^{ci +yi)) = .fi{x,yi) = 

= --x^ + -a;6 + -x"^ + yl- x^yi + x^yj + x^yi 

- 2x^yl + a;3y3 - x^i + ^xVi ~ x^vl + \^S\ + 2x^2/1 - ix^'yi 
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Figure 11: Newton Polygon for /i 

7' 




We call this new polynomial fi(x,yi), so now that we have our new function 
fi{x,yi), and our new goal is to find C2- We find the new Newton Polygon 
(Figure 11), from which we get /32 = 4 and 72 = 2. Similar to (3), we then have 



2/1 =a;^(c2 + 2/2) 



(7) 



where C2 is the coefficient of the second term in the series and 2/2 is 



j/2 = csx'^" 



C4X 



73+74 



(8) 



We substitute (7) into /i to get 

' = - ^.T^ + ^X' + \x^ + X\C2 + y2f - X\C2 + y2) + x\c2 + y2f 

+ a;^(c2 + 2/2) - 2a;^(c2 + ^2)^ + x^{c2 + y2f - x^{c2 + ya) 

+ -x^{c2 + y2f - xi°(c2 + y2f + ^a;i2(c2 + ya)^ + 2x^(c2 + 2/2) 

- 3x^(02 + y2)' + 2a;"(c2 + 2/2)' - \x^^{c2 + ^2)" - ^2;«(c2 + ^2) 

+ 5a;iO(c2 + y2f ~ bx^\c2 + ^2)' + \x^\c2 + ^2)' - \x^%C2 + y2)' (9) 



Again, lowest terms must cancel so, since [32 = 4, we see that | + C2 — C2 = 0. 
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Solving for C2 gives us C2 = |, and we substitute this C2 into (9) and divide by 
x~^ giving 



X ^f{x,x'^{c2+V2))^ h{x,V2) 



-x^ + -x® + -x'^ 



32 64 64 



^a;io + -a;2 - -a;^ + -x^ + y| + a;2y2 + a;^^ _ 2x^yl - -x^y2 
-x^yi + x^'yl - x^y2 + -a;^^^ + — a^^ys + 22;®y2 ^ x'^yl - — x^^ya 



„9„,2 ^9„,3 -'^^9„,4 I ^„10 
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,10„,2 



„10„,3 



„10„,4 



^x-y^2 - ^'Vi - 2^2/2 + 4^^'y2 + y.T^^y2^ + bx^'^yi + -x^'^y'^ 



5 12 5 ]^2 2 5 ]^2 3 5 12 4 1 12 5 

^x 2/2 - gS; ^2 - ^a; ^2 - ^2; y2 - -x 2/2 



(10) 



Figure 12: Newton Polygon for /2 



9 
8 

7 
6 

5 

4 
3 



f • • • 

• • • 
# • • 



K 



o 123 456789 10 U12 13 

We must now find C3. From the Newton Polygon for /2 (Figure 12), we can see 
that (33 — 2 and 73 = 1. Therefore, 



2/2 = x{c3 + ys) 



(11) 



where C3 is the coefficient of the third term in the series and y^ is 

y3 = C4a;T* + C5x'^''+'^= + • • • 



(12) 
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Substituting (11) into /2 we get 

15 57 

=^a:i2(c3 + 2/3)2 _ -a;i"(c3 + yg)' - 298a;9(c3 +j, 3) - x^\c:, + y^f 

- 2x5(c3 + y^f - 4^"(c3 + Vsf - 2;^(c3 + y:if + a:'(c3 + ys) 

+ 2^'(C3 + ys) - 42;i6(c3 + y3)4 _ _a;^(c3 + ^3)2 + -x\c^ + ^3)2 

- \x^\c, + ^3) + \x^\c, + y3)^ - y^%c, + y3) - ^^"(ca + y,f 

+ ^a;"(c3 + ya) + 3a;9(c3 + yg)' + ^^^^i^i + ^s)^ + a;"(c3 + ya)' 

3 15 

- a;^(c3 + ys) + 22;^(c3 + y?.f - 2^^^(^3 + 2^3)'"^ - 8^^''('=3 + ya)^ 

- 2^"(c3 + ya)' + -^x^ - g2;4 + -x3 - -x'^ic^ + y3) - -^x^^{cz + ys)' 
+ 2xi0(c3 + y^f - ^x^ + ^a;6 + ^x^ - ^x^ - ^x^ + ^x7(c3 + yg) 

- ^^12 + Aa;i° + x\c, + y,f + x2(c3 + y^f + 5x'Hcs + y^f (13) 

Since (3^ = 2 (which can be seen from Figure 12), when lowest terms cancel, 
we find that C3 + | + c§ = so C3 = — | . We substitute this value for C3 into 
(12) and divide by x^^ to get 

'h{x, a;(c3 + y3)) = Jz{x, y^) = 

1 ,„ 2 21 4 5 3 9 . 39 6 85 7 1 q 

- r - ''^ + y ^ - r " r ~ u^ + 32^ - r 

+ fx'Vs f -^°y3 + fx'Vs - y x«y3 + f x^^ya - l^Vs + QxVs 
" 2^^3 + ya;«ya - yx^y3 + -x^y32 - — x*^y32 + -x^% - -^^^V^ 
+ ^x"yf - Ix'V^ + ^-"ya - 5x^2^33 + ^.12^4 ^ 2^Syi _ 5^6^^ 
+ 2x6y2 - ^x"y32 + 6x"yf ~ ^x'"yl - ^x% + ^x^y^ - x2y3 + ^2^32 

- ^2;"y3 + —x^^yl - -x"y3 + x^y^ - 2x3y| - x^y^ - -x"y^ 
+ \x^%l + x^yl - y^^y, + Ix'^yl - \x'^yl + ^x^Ys ^ ^^^^yf (14) 



we 



X-2 
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Figure 13: Newton Polygon for f^ 



o 1 23 45 67 89 10 U12 13 14 15 16 



To find C4, we repeat the process for f^. Looking at the Newton Polygon (Figure 
13), we see that ^4 = 1 and 74 = |. From this we have 



2/3 = a;2(c4 + y4) 



(15) 



where C4 is the coefficient of the fourth term in the Puiseux series and 2/4 is 



y4 = C^x^^ + C5X''''+^'- + 



(16) 
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Substituting (15) into /a, we get 

,^ 3 c, , ,9 123 23, , 1 _ -tn, ,9 11 19/ N 

^-x^Ci + ViY + -—X 2 (c4 + Vi) + -x + 6x'"(c4 + 2/4) - ^x 2 (c4 + 2;4) 
z oz Z o 

15 27 N „ 4/ n9 „ 19/ n'^ 5 31 , , 

- TE^ ' (c4 + Vi) - 2a;*(c4 + y^Y + 2x 2 (04 + y^y - —x ^ (04 + 2/4) 

16 32 

- -x¥(c4 + 2/4)^ + a;5(c4 + j/4) - a:;¥(c4 + ^4)^ + -a;i^(c4 + 2/4)^ 

5 29 N 9 8/ n9 1 T^/ n4 1 35/ s r; 

+ gs; 2 (c4 + 2/4) + -x«(c4 + y^Y - i^x'^c^ + Vi) - ^^ ^ (c4 + ^4)" 
+ ya;"(c4 + y4)^- -a;^(c4 + y4) + -^x^'^^c^^ y^f - hx^ {c^ + y^"^ 

59 15 15 9 81 n, ,9 -/ ^ 

- — a; 2 (c4 + j/4) - — x2 (c4 + y4) - —x (C4 + ^4) - x^ (04 + ^4) 

- 2^^ ' (C4 + V^f + 22^(^4 + y^T - -jx^ {ci + yiY - -^x 2 (c4 + yYY 
+ —XT (c4 + ^4) + yxl3(c4 + y4)2 + x^ {c^ + ^4)^ - -x^(c4 + y^Y 

+ y xl4(c4 + yYY - -x2 + — x^ ^ -x3 + — X^ (04 + yY) 

9 21 o 5io/ ,4 9c, 39 c 85 7 5i7, ,4 

- 2^^ (^4 + VY - -^^ {Ci + y4)* - -^X'> - — X^ + — X^ + -X1^(C4 + 2/4)^ 

37 21 . , In 7 o 5 31 . , q 17 in 

xT(c4 + W4) X-* X** + -XT(c4 + V4)'^ H X^" 

8 ^ 4 -r y4; ^ ^ -r 2 V 4 -r y4; -r ^^ 

x" H x^^ x"h x" + 2x'^(c4 + ?y4)^ x^(c4 + yY 

64 64 64 32 ^ '^^ 8 ^ ^ -r y4; 

x3(c4 + 2/4)^ + 6x¥(c4 + yY^ - -^x^^a + yY^ + x(c4 + 2/4)^ (17) 

8 



-x" 



Since /34 = 1, when lowest terms cancel, we find that 5+C4 = so C4 = ± 2^- 
So the two solutions given by the Newton Puiseux Algorithm are 



y = -x^ 



^x2 



x^ x^ (-2x)^ 

y " y ^ 64~ 



x^ (-2x)^ 



2 2 64 



as one of the jets. It can be shown that the the possible values of C4 are equal 

11 

to the roots of ^^-g^^ — , and our solutions for C4 seem to be imaginary. However, 
x^ forms a cusp, and the fact that there is a — x in the jet simply means that 
the cusp opens left, so when x is negative, the coefficient of the -^ term in the 
jet is real. 
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